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NEW POLYTOPE DECOMPOSITIONS AND EULER-MACLAURIN 
FORMULAS FOR SIMPLE INTEGRAL POLYTOPES 

JOSE AGAPITO AND LEONOR GODINHO 


Abstract. We use a version of localization in equivariant cohomology for the 
norm-square of the moment map, described by Paradan, to give several weighted 
decompositions for simple polytopes. As an application, we study Euler-Maclaurin 
formulas. 


1. Introduction 

The interplay between symplectic geometry and combinatorics through the study 
of moment maps and the use of equivariant cohomology and the geometry of toric 
varieties is a well known and fertile theme in mathematics. See for example the 
work of Brion-Vergne Fvn . Cappell-Shaneson iran . fcsi . Ginzburg- Guillemin- 
Karshon |GGKj. Guillemin |Gul|. IGu2l. Morelli !H1, and Pommershein and Thomas 

eu, m 

In this paper we use a version of localization in equivariant cohomology for the 
norm-square of the moment map, due to Paradan (see 0), to motivate several new 
weighted polytope decomposition formulas. Indeed, in Section |21 applying this local¬ 
ization principle to the particular case of toric manifolds, we obtain weighted polytope 
decompositions for Delzant polytopes |Del (see Example IT7I) and then, by a purely 
combinatorial argument, we show, in Section 0] that they are in fact valid for any 
simple * 1 polytope not necessarily the moment map image of a toric manifold (cf. The¬ 
orem mi Moreover, still in Section QJ we use these decompositions to obtain new 
ones iTheorem l4.2ll that generalize both the Lawrence-Varchenko decomposition (see 
m and Cj) and the Brianchon-Gram formula (see E3, E3, E3 and |Grj l (cf. Re¬ 
mark 138ft . 

The well-known classical polytope decomposition formula of Brianchon-Gram ex¬ 
presses the characteristic function lp of a convex polytope P as the alternating sum 
of the characteristic functions of all tangent cones to the faces of P. By flipping the 
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edge vectors emanating from each vertex of P in a systematic way using a polarizing 
vector, we obtain the Lawrence-Varchenko decomposition (also known as polar de¬ 
composition) which expresses the characteristic function of a convex simple polytope 
(only) in terms of the characteristic functions of polarized cones supported at the 
vertices. Karslion. Sternbere' and Weitsman IKSWI and Aeu/nito EU gave weighted 
versions of this decomposition by assigning weights to the faces of the polytope and 
of the cones in a consistent way. Our polytope decompositions combine the above two 
formulas. Like Brianchon-Gram they express lp in terms of characteristic functions 
of cones with apex the different faces of the polytope. However, these cones may no 
longer be the ordinary tangent cones to the polytope. Indeed, in our formulas, to 
each face is assigned a different polarizing vector and we flip the edges of the tangent 
cones accordingly. In the first decomposition formula (Theorem 14.11) these polarizing 
vectors are obtained by choosing a suitable starting point e (the same for all vectors) 
and then taking as end points its orthogonal projections /3(e, F) onto the faces F of 
the polytope, whenever these projections are nonempty (cf. Figure [2J|. In the second 
decomposition formula ("Theorem 14.211 we take the vectors e — /3(e,F) as polarizing 
vectors (instead of (3(e,F) — e as above). This second formula generalizes both the 
Lawrence-Varchenko and the Brianchon-Gram relations in the sense that, choosing e 
in suitable regions of we obtain these polytope decompositions. 

As an application, in Section[5l we use the new decompositions to give new Euler- 
Maclaurin formulas with remainder similar to those of Karshon-Sternberg-Weitsman 
|KSWj and Agapito-Weitsman I3WI . The classical Euler-Maclaurin formula computes 
the sum of the values of a function / over the integer points of an interval in terms 
of the integral of / over variations of that interval. This formula was generalized by 
Khovanskii and Pukhlikov (see EE2 and |KP2l l to a formula for the sum of the values 
of an exponential or polynomial function on the lattice points of a regular integral 
polytope, by Cappell and Shaneson IGSTI . IGS21 . 0, Guillemin lGn2l and Brion 
and Vergne IBrVlI . to simple integral polytopes, and by Berline, Brion, Szenes and 
Vergne nwsi . UT7v2i Eg, to any rational polytope. Note that all these formulas are 
exact and valid for sums of values of exponential or polynomial functions. Moreover, 
the formula in [Bill has the additional feature that it is local , in the sense that it 
is given as a sum of integrals over the faces of the polytope of maps D(F ) ■ p, for 
operators D(F) depending only on a neighborhood of a generic point of the face F. 
In IKSWI . Karshon, Sternberg and Weitsman prove a formula with remainder for the 
sum of an arbitrary smooth function / of compact support, on the integer points of 
a simple polytope. There, the remainder is given as a sum over the vertices of the 
polytope, of integrals over cones with those vertices, of bounded periodic functions 
times several partial derivatives of /. In our formulation, both the Euler-Maclaurin 
formula and the remainder are given as a sum over all faces of the polytope (not only 
over vertices) of integrals over cones with apex the affine spaces generated by those 
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faces (see Theorem 15. 111 . Moreover, since our formula generalizes to symbols (in the 
sense of Hormander jHj), we show that, in the case of polynomial functions, we also 
obtain an exact Euler-Maclaurin formula for the sum of the values of a polynomial 
function p over the integer points of the polytope. This relation is a weighted version 
of the exact Euler Maclaurin formulas obtained in |BrVT| and in |KSW| . 

2. Critical set of the function ||/z £ || 2 

Let ( M,u >) be a compact connected symplectic manifold equipped with a Hamil¬ 
tonian action of a torus T. Denoting by t the Lie algebra of T and by t* its dual 
space, we consider the moment map fi : M —> t* associated to this action (that is, the 
T-equivariant map determined, up to a constant, by the equation d(p,X) = l(Xm)cj 
for all X £ t). Since the action of T on t is trivial, the perturbed map p e := fi — e for 
e e t* is also a moment map for the action of T. 

Hereafter we will choose a scalar product on t* (which induces a linear isomor¬ 
phism j : t —> t*) and consider two kinds of orthogonality on t*: the orthogonality 
resulting from the duality between t and t* and the j -orthogonality defined by the 
scalar product. We will begin by reviewing the structure of the critical points of the 
moment map and then define an index set B which will enable us to define a partition 
of Cr(||/i £ || 2 ), the set of critical points of the function ||/z £ || 2 . 

Following Paradan in |P|, we will consider a slightly modified definition of critical 
point. By the usual definition, a point a; is a critical point of the moment map /i if 
its stabilizer, Stab(x), contains a subtorus of dimension 1. However, if the action is 
not effective, all points of M will be critical points of fi. To avoid this situation we 
take the subgroup Sm '■= n x eM Stab(a;), called the generic stabilizer, and make the 
following definition: 

Definition 1. The critical points of the moment map fi : M — > t* are the points 
x € M for which Stab(a;)/iS , M is not finite. 

Let T' be a subtorus of T containing Sm such that T'/Sm is not finite. Every 
connected component F' of M T is a symplectic submanifold of M and P := p(F') 
is a convex polytope in t* equal to the convex hull of the image of the fixed points of 
T contained in F' (c.f. m and ESI). Moreover, the Lie algebra of T' is included 
in P^ (the set of vectors of t orthogonal to P), and, denoting by Tp the subtorus of 
T generated by exp(P J -), we have that the Lie algebra tp of Tp is equal to P- 1 . The 
manifold F' is then a connected component of M Tp where T/Tp acts quasi-freely, 
that is, the generic stabilizer of F' is a subgroup of T and its identity component is 
equal to Tp. Knowing this and denoting by Aff(P) the affine subspace of t* generated 
by P, we consider the following sets: 

(2) B' := {convex polytopes P C t* for which there exists a connected 
component F' of M Tp with fi(F') = P}~ 
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(3) 

B := 

= (Aff(P) I P £ B'}- 



(4) 

B’a'. 

= {P £ B' 1 P C A 

and dirn(P) < dim(A)}, for A £ B; 

(5) 

Cg 

: = t* \ u P, 

^reg 

:=A\ (J P; 



P€B'\h(M) 


pgb^ 

(6) 


Wa A reg + j(tA) ; 

for A £ B; 

(7) 


w = n 

Wa ■ 



AeB 


Note that the set B' contains all the faces of the polytope fi(M) and that, if P is a 
polytope in B ', all its faces are also in B'. Again, for A £ B, we will denote by T/\ 
the subtorus of T generated by exp(A x ). With this notation we have the following 
proposition which characterizes the critical set of ||/4 e || 2 (c.f. jK] and [P] for details): 

Proposition 8. (Kirwan) For every e £ t*, 

Cr(||Me|| 2 )= |J M Ta (~l ^i _1 (/?(e, A)), 

AeB 

where, for an affine subspace A oft*, /3(e, A) is the orthogonal projection of s on A. 
Moreover, for every e £ Wa, the set 

(9) Ca := M Ta n /z -1 (/?(e, A)) 

is a submanifold of Ad on which T/Ta acts locally freely, the set W is dense in t* and, 
for every £ £ W, the submanifolds C\, for A £ B, form a partition o/Cr(||/z e || 2 ). 

Since the group T/Ta acts locally freely on the manifold C%, we can define the 
quotient := C^/(T/Ta) which will be an orbifold. Moreover, for each connected 
component F of C^, we will consider the subgroup S A (F) := C\ xe p Stab(x) and the 
map F i—> |S' A (i ;l )| (defining a locally constant function on M ^) which will be called 
\S A \. 

Remark 10. Note that, if dim(T) = l/2dim(M) (i.e. if M is a toric manifold), then 
M = ^ _1 (/3(e, A ))/T is either empty or a point. 

From now on we will consider e £ W. 

3. Localization Formulas 

3.1. Equivariant cohomology. Let us begin by reviewing the different T-equivariant 
de Rham complexes on M. We have three spaces of equivariant differential forms on 
M, f Vf{M) C fl“(M) C fly°°(M), respectively with polynomial, smooth and gen¬ 
eralized coefficients. The model fi^(M) is due to Cartan while the other two were 
studied by Berline, Duflo, Kumar and Vergne (see IBP V| . rBV2j . |BVj and Kvj ) . On 
f2“(M) we have the differential dx defined by 

(d T a)(X) := (d-i(X M )){a(X)) 





NEW POLYTOPE DECOMPOSITIONS AND EULER-MACLAURIN FORMULAS 


5 


for every a G and X G t. The corresponding de Rham T-equivariant com¬ 
plexes on M are defined as := (fl^? (M),dr) and H^(M) := (f l^(M),dr), 

and called the T-equivariant cohomology with C°° and polynomial coefficients. More¬ 
over, the differential dr defined on extends to Q^°°(M). Indeed, if {e 1 ,..., e r } 

is a basis for t then, for every r) G f2!j?(M), we have 

r 

{d T {v),<t >) ~d{v,<t>) ~ X k ), 

k= 1 

for every compactly supported density </> on t, where X \,..., X r are the coordi¬ 
nate functions on t. The corresponding de Rham T-equivariant complex is denoted 
:= dr) and called the T-equivariant cohomology on M with 

generalized coefficients. We can also consider equivariant classes of compact sup¬ 
port obtaining T-equivariant compact support de Rham complexes: t (M) := 

(I4 iCpt (M),d T ), H^ cpt (M) := (fi§? cpt (M),dr) and := , d T ). 

3.2. Equivariant Euler classes. Let p : E —► M be an oriented T-bundle and let 
p* : fl cpt (E) —» f 1{M) be integration along the fibers. There is a unique equivariant 
class u G Ht cpt(E) such that p*w = 1 on M called the equivariant Thom class of E 
and denoted by Thom'r(T). Its restriction to M is the equivariant Euler class of the 
bundle E, i.e. ex{E) := i*ThoniT{E), where i : M —► E is the inclusion map. 

Let us assume now that there is an element f3 G t for which the zero set of the 
vector field on E generated by (3 is equal to M and let Tp be its stabilizer in T. 
Then there is a T^-equivariant class with generalized coefficients e^ 1 (E) such that 
e^ 1 (E)eT(E) = 1 (see fP( for details). 

Example 11. Consider the trivial bundle E := M x C equipped with a T-action 
which is trivial on M and which, on C, is determined by the weight a G t* (that is, 
exp(X) ■ z := e l ^, x ) z) for X G t and z G C). Choosing /3 G t such that (a,/3) ^ 0, 
we have er(E){X) = —^(a,X) and, “polarizing”, that is, taking a + := ep a with 
(a + ,/3) > 0 and ep = ±1, we have 

POO 

e~ 3 1 {E)(X) = 2mep e i < a+ ' jr > * dt 
J o 

as generalized functions. Taking the Fourier transform we obtain the equality of mea¬ 
sures on i*. (E)) = 27tz ep H a +, where H a + is the Heaviside measure associated 

to a + defined by 

POO 

(H a +,cj)) = / 4>(ua + )du, 

Jo 

for every <fi in the Schwartz space of rapidly decreasing functions on M. 

Example 12. If M = {F} is a single point, fixed by the action of T, the bundle E 
decomposes as a sum of non-trivial 2-dimensional real representations of T, E := 
L\ © • • • ® Lk —► F, with the action of T on each Lj determined by a weight ay G t*. 
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Following Paradan (cf. Proposition 4.8 in |P]) we obtain the expression for the Fourier 
transform of e(j 1 (E): 

F( e p 1 (E)) = (2ni) k e 0 H at *---*H at , 

where we polarize each a 0 according to some P £ t (such that (a,, (3) yt 0 for j = 
1,..., k), obtaining := ay with = ±1, and we take ep := ] _ [ / ( = i e^. Note that 
* denotes the convolution product. This measure, supported on the cone K + af + 
-h M+a^, is defined by 

n OO 1*00 fe 

{H+ *■•■* H a + ,(/>) = / ■■•/ (t>{S^Uial)dui...d,Uki 

1 k Jo Jo 

for every rapidly decreasing function on M. 

3.3. Localization. Using the sets B , W and C A defined in (JSJl, Q and © of Sec- 
tionO and the orbifold M A = C A / (T/Ta), Paradan proves the following localization 
theorem: 


Theorem 3.1. (Paradan) Let e £ W and let rj £ Qj?(M) be a closed form, 
on C~°°( t) we have 


M 


= E 


Aee 


where I A (rj) is the generalized function supported on tA defined by 


Then, 


(13) 


where 


r A (rj)(X 1+ X 2 ) = 

= (2™) dimA f J ^ K -k A (r 1 )(X 1 )e^(E A )(X 1 )od(X 2 -w A ), 


(i) the variables Xi,X 2 are respectively in tA and i/i A (note that, for each A £ B, 
t decomposes as a sum of vector spaces t A and t/t A , where i A is the Lie algebra 
of the subtorus T A generated by exp(A J -)j; 

(ii) k A : Hf?(M) — > Hf?^(M A ) is the Kirwan map (see lK]); 

(hi) Pa '■= j _1 (/3(e, A) — e), where P(e, A) is the orthogonal projection of £ on A; 

(iv) E a := N a /(T/T a ), where N A is the normal bundle of M Ta inside M, re¬ 
stricted to C A ; 

(v) the operator o: 


^(CDxn-jci) - ^°°(C7i) 

(r),i') i—> rj o v 


is defined by 

(V o v, <j>(X)dX) := (r), (v, <j)(X i + X 2 )dXj)dX 2 ) 

for every density 4>(X) of compact support on t; 

(vi) w A is the equivariant curvature of the principal orbibundle C A —> M A ; 
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(vii) the equivariant form S(X 2 — w A ) £ f2 T ?^(M^) is defined by 

(8(X 2 - w A ),(f>{X 2 )dX 2 ) = <f>(w A )vol(T/T A ,dX 2 ), 

for every function 4> £ C°°( t/tA), where vo1(T/Ta, dX 2 ) is the volume of the 
group T/T A with respect to the Haar measure compatible with dX 2 . 


Example 14. If A = {p} is a vertex of the polytope p{M), then C A = p, l (p) is a 
connected component F of M T and 

I\ p} mX) = f i* F (v)(X)e^(N F )(X), 

J F 

where N F is the normal bundle of F inside M and /3 P = j l {p — £)• If, in addition, 
the action of T is toric, then F is an isolated point. Moreover, taking ?/ = , where 

w* 1 is the equivariant symplectic form on M defined by u^(X) := u>(X) + (n,X), we 
obtain 


(15) 


F {p} (e i “’)(X) = e p e 


dim(M)/2 

i{P,X) 

j=l 




where the a ’s are the polarized weights of the action of T on the normal bundle of 
F (i.e. the polarized edge vectors at p) and e p is the sign obtained by polarization. 
Taking its Fourier transform we get 

(16) HI\ P} (e“ # )) =e p 8 p *H+ i *-..*H+ n , 


where 6 P is the Dirac measure on p £ t*. Moreover, the measure m is supported on 
the polarized cone := p + M + af + • • ■ + R + a+. 


Example 17. For a toric manifold M 2n with moment map pL and for ?/ = where 
again is the equivariant symplectic form on M (cf. Example [EJ, the reduced space 
M e A for e € t* is either empty or a single point. Moreover, denoting by /3i(e, A) the 
orthogonal projection of /3(e, A) onto t^, formula ltl3l) becomes 


(18) 


I e A (X 1+ X 2 ) 


= (2m) 6 a 


s i(/3i(e,A),Xi) 

W\ 



o S 0 (X 2 ), 


whenever /3(e, A) fl /i(M) is nonempty, where e A = TIy=i e A with a\ • = a A ,j is 
the sign obtained by polarization, where r A is the codimension of A, where IS 1 ^ is 
the order of the orbifold structure group of the point in M ^ inside the toric orbifold 
E a := N A /T/T Al and where the -’s are the polarized weights of the action of T A 
on E a restricted to the normal orbibundle of the fixed point. Equivalently, |5 A | is 
the order of the orbifold structure group of the point in inside the reduced space 
I t T/T A (^ 2 ( £ > A))/T/Ta, where ^t/t a is the moment map for the T/TA-action on M, 
and /3 2 (e,A) is the orthogonal projection of /3(e,A) onto (t/tA)*- 
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Taking the Fourier transform of © we obtain 
(19) ( 27r *) n |^| ((fyi(e,A) * H a +^ *■■■* H ai J o l (t /u)*) , 

which is supported on the polarized cone C^ £ := (3i(e, A) + R + a A 1 + ■ • • + 

R + a A . Moreover, changing variables, we obtain 


( H a + , * • • • * H a+ > </>) = 


/•OO r A 

/ </>(S2 Ui a^i) dui ... du rA = 

do 


I det (aXi)i\Jc* |det(at i ) i | 0 


5 4*) 5 


for any rapidly decreasing function </>, where is the characteristic function of the 

'■-'o 

cone Cq := R + a A 1 H-bR + a A . However, since the oa/s are the weights of the 

action of Xa on the toric orbifold E& at the fixed point, we have | det (a A | = 
and so C3 becomes (27n) ra eA l c tt o l(t/t A )*- Indeed, denoting by £* the lattice dual 
to the weight lattice l of Xa, the orbifold structure group F of the fixed point is 
isomorphic to £/£, where £ is the lattice of circle subgroups of Xa (cf. m for details), 
implying that 


l*5 ,A | = |F| = | det X| 


1 _ 1 

I det i\ | det (aJ^),| 


(see for example El)- On the other hand, since the Fourier transform of f M e tuJ,t is 
the direct image /i,(^ ro i) of the Liouville measure dniL '■= ui n /n\ on M (which is 
supported on /x(M)), we obtain (modulo (2ni) n ) 


(20) 1 ii(M) = A ) eA 1 ci ’ 

A SB 

where ^>(e, A) is equal to 1 when /3(e, A) 0 /x(M) is nonempty, and zero otherwise. 
Note that Formula is valid up to boundary effects on the polytope /x(M). 


4. Polytope decompositions 

In this section wc will show that the polytope decomposition for Delzant polytopes 
that was obtained in CT . remains valid for any compact convex simple polytope. 
Moreover, we will give a weighted version of this decomposition that also holds on 
the boundary of the polytope. 

Hereafter, we will consider the usual Euclidean inner product (,) of Let P be 
a compact convex simple polytope in R d and let B' be the set of faces of P. For each 
F £ B 1 we write Ap for the affine subspace of R. d generated by F. Then, just as in 
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Section |3 we have the following sets: 

(21) B := {A f | F e B'}\ 

(22) B' a := {F e S' | F C A and dim(F) < dim(A)}, for A G 6; 

(23) A reg := A \ U F: 

(24) Wa := A reg + A 1 - for A 6 B; 

(25) W:= p| Vh A - 

AeB 

The set IT is a disjoint union of open sets which we will call Paradan regions (see 
Figure |3 for an illustration). In fact, W is the complement in of a finite set of 
walls of codimension 1, W c = E\ U • • • U Ek U {facets of P}, where each wall Ei is 
contained in a hyperplane perpendicular to a family of elements of B. Note that, in 
the case of a moment polytope /x(M), the set W is the same as in 0. 



Figure 1. Paradan regions for a triangle. 


4.1. Tangent cones. For each A £ B, we define the tangent cone of P at A by 

C A := {y + r(x - y) \ r > 0, y £ A, x G P}. 

It is a full dimensional cone with apex A (i.e. A is the maximal affine space contained 
in Ca). Taking e in R d , we denote by /3(e, A) its orthogonal projection onto the affine 
subspace A, and we take the intersection of /3(e, A) + A- 1 (the orthogonal space of A 
at (3(e, A)) with the tangent cone of A, 

C A y £ :=#,A)+A i )nC A , 
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which is now a pointed cone 2 with vertex (3(e, A). The cone Ca is the direct product 
of the affine space A and the pointed cone C A x e . Then, considering vectors £ 
R d along the edges of C a _l )E , pointing away from the vertex ( j = 1,... ,r A , where 
r A = dimC A x £ = codim A), the tangent cone C A can be written as 

ra 

(26) C A = A + C A r E = A + 'y ' R+CkAj, 

3 =1 

that is, C A is the cone along A which contains P and is bounded by the affine spaces 
in B which contain A. The vectors a A y (which are only determined up to a positive 
scalar) will be called the generators of C A (see Figure |21 for an illustration). 


/ ..••• \ 


l • ' ’ ' ^ jS*' I 


(3{e, &ab) * 


A 

P(", 


.0 

"A" 


*»/3(e, Aac) 


I , v & 


B 




Figure 2. Projections and generating vectors for some faces of a triangle. 


4.2. Polarization. We will “polarize” the tangent cones in the following way: let 
us consider a point e £ W; as a direct consequence of the definition of W, for each 
A £ B we have (/3 A . a A y) 7 ^ 0 for all j = 1,..., r A , where /3 A := /3(e, A)— e; polarizing 
the vectors a Aj according to /3 A , that is, taking the vectors a\ j = e^ A a Aj with 
(a\ p /3 a) > 0 and e^ A = ±1, we define C^, the polarized tangent cone of P at A by 

r a 

(27) C“,:=A + ]Tr+ Q + .. 

3 =1 


2 


A cone with a single point as apex (also called vertex). 
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4.3. Weighted characteristic functions. Let us now see how to assign weights to 
each affine space in B in order to obtain a weighted version of ®) that also holds for 
points in dP. Take Ai,..., A n, the codimension-1 elements of B , that is, the affine 
subspaces generated by the facets of P. Each A E 8 that is generated by a non-trivial 
face of P (i.e. such that A ^ 0 ,K. d ) can be described as an intersection 


n 


where Ja denotes the index set of the hyperplanes A, that contain A. Note that, 
since P is simple, the number of elements of Ja is equal to ta, the codimension of 
A. To each A,; we assign an arbitrary complex number q and to each affine space 
A e B\U. d we assign the value ( h' Moreover, to A = we assign the value 1. 

This amounts to defining a weighted function w : R d —> C, given by w(x) — riiej Air <7i) 
where A x is the smallest-dimension element of B that contains x. With this function, 
we define the weighted characteristic function 


(28) 


1 p(x) 


w(x), if x e P 
0 , otherwise 


Remark 29. Similarly, for each A £ £>, we define weighted characteristic functions for 
the tangent cone Ca and for the polarized tangent cone Cjj^, but now assigning the 
weight 1 — qi to each facet of Ca that is flipped in the polarization process. 


The fact that this assignment of weights fits together and makes our polytope 
decomposition hold even for points in dP, relies heavily on the following Lemma as 
we will see later in the proof of Theorem im 

Lemma 30. Let I = {1,..., n}. For any qi £ C with i £ I, we have 

(si) n %=i+ 55 life 

iei 0 ^Jci jeJ 

Proof. We will use induction on the cardinality n := (I) of I. Clearly, formula (13111 

is trivial for n = 1. Let us then assume that ED holds for the set I' = {1, ..., n — 1} 
and show that it holds for I = {1 ,,n}. Indeed we have 

n qt = (ii qt ) qn =( 1 + 55 n fe _ ^) fe ~ 1 +^ 

iei \t=1 / \ e^J'cr jeJ 1 ) 

= 1 + {in - 1) + 55 ( II fe ~ ^)fe “ 1) + IT fe - 1) ) ■ 

0/j'c/' \jeJ' ieJ' ) 

However, on the other hand, 

life -!)= 55 (9n-1) n fe- i )+ 55 nfe _i )> 

0/jc/ jeJ Jci j£J\{n} 0 ^Jci' ieJ 

with nGJ 


and the result follows. 


□ 
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4.4. Decomposition formulas. Finally, defining <p(e, A) as 


(32) 


¥>(£> A ) 


1 if /3(e, A)flP /0 
0 otherwise 


we obtain, for each e £ W, the following polytope decomposition formula (See FigureEl 
for an illustration). Note that, by the definition of ip, this formula only takes into 
account the polarized cones for which (3(e, A) fl P / 0 , that is, those for which 
the orthogonal projection of e onto A is in P. 


9293 



9193 9l 9192 


9293 



Decomposition corresponding to e in region 7^2 of Figure 2. 



Decomposition corresponding to e in region TZi of Figure 2. 


Figure 3. Polytope decomposition for a triangle 


Theorem 4.1. For any compact convex simple polytope P of dimension d in and 
for any £ in W, we have 

(33) l?=E(-ir^( e »A)lS>. 

A £B A 

where the sum is taken over the set B of affine spaces generated by the faces of P, 
where is the tangent cone of P at A £ B polarized with respect to the vector 
/3(e, A) — e (where /3(s, A) is the orthogonal projection of £ onto A), where m/\ is 
the number of generators of the tangent cone Ca whose signs change by polarization, 
and where lp and 1“ # are the weighted characteristic functions of the polytope and 
of the polarized cones respectively. 
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Proof. We will prove this formula in two steps. First, for each x we will find an £ in a 
Paradan region for which formula El holds. Then, we will show that the right hand 
side is independent of the choice of e. 

Step 1: Suppose x P. If x is in W then we choose e on the same Paradan region 
as x. If, on the other hand, x is in the complement of W, we choose e on any Paradan 
region contiguous to x inside the complement of P. In both cases, none of the cones 
on the right hand side of 13311 contains x and we obtain 

0 = 1?W = E (-irMe, A)l“, (x) = 0. 

AeB 

If x is on the boundary of P, we choose e in the interior of P close enough to x so 
that ip(e, A) = 1 for every A e B that contains x. In this case, all polarized cones 
that contain x are pointing away from e. Moreover, all generators of the corresponding 
tangent cones Ca are flipped and so to a = codim(A). Denoting by A x the smallest 
dimensional affine subspace in B that contains x, we have lp(x) = Y\,^j& while 
the right hand side of El is 

RHS = £ (—ir>( £ ,A)l”,W= £ (-l) m *w A (x) 

AeB x A As B x 

= i + £ (-i) mA n( 1 -®) = i + E II(®- 1 ) = 1 +E life- 1 )’ 

AeB x \R d A6B x \R d iGJA JC Ja x j€ J 

where B x C B is the subset of elements of B that contain x, where u> A (x) := rw 1 
q,j). and where we used the fact that to a = codimA = #J A and that A £ B x iff 
J A C Ja x - Consequently, from Lemma 1301 we conclude that 

i?(x) = E(-irM e >A)is<(x)- 

AeB A 

If x is in the interior of P then, taking e = x, none of the polarized in 1331 
contains x except the one corresponding to the top dimensional face P of the polytope. 
In this case, </?(£,R d ) = 1, and, formula ESI evaluated at x gives 

l = l£(x) = ^( e ,IR d )l£ # (*) = !• 

R d 

Step 2 : Recall that the complement W c of IT is a finite family of walls of codimension 
1, W c = Ei U • • • U Ek U {facets of P}, where each wall Ei is perpendicular to a family 
Bi of affine spaces contained in B. Let eq and £2 be in two contiguous Paradan regions 
TZi and IZ 2 respectively, and let E be its common “wall” (either one of the Efs or a 
facet of P). Let e* be any path in from £1 to £2 that crosses a single wall (i.e. E) 
once. When e t crosses E, the sign of 

(34) (P(e t , A) -£t,a A) fc) 

(for A £ B \ and for a generator a A: fc of the tangent cone at A) flips exactly when 
AncfP is contained in E and aA,k is perpendicular to E. Hence, if dim A ^ d— 1, this 
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cube, with e\ and £2 in two contiguous Paradan regions. 


sign flips iff A is contained in (exactly) one affine space A of B % with dim A = dim A+l 
(see Figure 0] for an illustration). Indeed, we can take 

A := A + spanaA,fc & Bi 

and unicity follows from dimensional reasons. On the other hand, if dim A = d — 1, 
the sign of m flips iff A fl P C E. In this case, we define A to be the entire space 
R d . 

Let us assume without loss of generality that the sign of m flips from negative 
to positive as e t crosses E. In this case, the polarized tangent cones at A before and 
after e t crosses the wall are 

(C$0 1 = A + ^M+a^.-R+a A ,fc and (C$J 2 = A + ^R+a^ + R + a A , fc . 

j^k j^k 

Hence, the corresponding contributions of A to the right hand side of El are 


±1 


w 



Tl 


(exp¬ 


and 
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Figure 5. Polarized tangent cones at the affine subspaces Ajj, 
A dh, Aim and A aehd for the cube in Figure 0] 


Note that the union of the two cones (Cj ^) 1 and (C^J 2 is the polarized tangent cone 
at A, C~ , for both e x and e 2 (cf. Figure 0, and so 

\ w , += r, 

( C ^) 1 + ( C ^) 2 1 c ^- 

On the other hand, we have A) OP/0, while j3(e 2, A) OP = 0. Hence, the 

corresponding contributions of A to the right hand side of 01 are 

=Fl^ti and 0 . 

Indeed, f3{E, A) n dP = An dP and /3(/3(ei, A), A) — A) = uaA,k for i = 1,2, 
with n < 0 and r 2 > 0 (cf. Figures 0 and 0. 

Consequently, the differences of the contributions of A to the formula in 01 before 
and after e t crosses the wall, and those of A, sum to zero. 

Moreover, for a given A £ B , if (p(et. A) changes when crossing E, the intersection 
of A with E contains A n P for (exactly) one element A of B with dim A = dim A — 1 
and the result follows. □ 

Remark 35. This new polytope decompositions 01 generalize the weighted version 
of the Lawrence-Varchenko relation for a simple polytope presented in HD There, 
the edge vectors emanating from each vertex are flipped in a systematic way using 
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a polarizing vector, and the weighted characteristic function of the polytope is ex¬ 
pressed (only) in terms of the weighted characteristic functions of the polarized cones 
supported at the vertices. In iTiali . not only the polarization is carried out differently, 
but, for some values of e, we consider the weighted characteristic functions of polar¬ 
ized tangent cones to faces other than vertices. Indeed, given e £ W, we obtain a 
different polarizing vector for each face of the polytope by taking e as starting point, 
and its projections onto the faces of the polytope as end points, whenever these pro¬ 
jections are nonempty. Then we polarize the tangent cones of the corresponding faces 
accordingly. 

4.5. Other decomposition formulas. If we polarize the generators of tangent cones 
with respect to e — /3(e, A) instead of /3(e, A) — e, and multiply each term on the right 
hand side of El by a factor (—l) dlmA , we obtain new polytope decompositions, 
under the same hypotheses and statements of Theorem rrn 

Theorem 4.2. For every compact convex simple polytope P of dimension d in W* 
and for any e £ W, we have 

(36) lp=J2 (-l) mA+dim Va A)l^, 

where the sum is taken over the set B of affine spaces generated by the faces of P, 

where is the polarized tangent cone of P at A £ B with respect to the vector 

e — /3(e, A) (where /3(e, A) is the orthogonal projection of e onto A), where m/\ is the 

number of generators of the cone Ca whose sign changes by polarization, and where 

1 f, and are the weighted characteristic functions of the polytope P and of the 
c A 

polarized cones respectively. 

Proof. The fact that the right hand side of does not depend on e can be proved 
as in the proof of Theorem rm Hence, we just have to show that we can find an e in 
some Paradan region for which EH holds. 

For that, let us choose an e such that <p(e, A) = 0 for every A £ B with dim A > 0, 
i.e. we choose 

e £ P| (F + A p) . 

F a facet of P 

Then formula E3> becomes 

(37) 1?= £ (-i)"*” 1 c*- 

v a vertex of P 

Choosing a vector £ £ R d such that, for every vertex v of P , > 0 whenever 

(e — v,a v j) > 0, where the ot v ,fs are the edge vectors at v, formula EU becomes 
a weighted version of the Lawrence-Varchenko polytope decomposition (see n ei 
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IkSWI and HU), where the tangent cones at vertices are polarized according to £, 
and the result follows 3 . 

This choice of polarizing vector £ can be done in the following way: first we consider 
the vertex vo of P that is furthest away from e; clearly, for vo we have C^ 0 = Cjj, 
(where this cone is polarized with respect to the vector e — vo)' : then, for any other 
vertex v and for each edge vector a v j satisfying (e — v,a v j) > 0 , we take the hy¬ 
perplane H 3 • through vo which is perpendicular to a v j ; these hyperplanes intersect 
at vo and each of them separates the whole space R d into two open regions. Let us 
denote by (iJj ■)+ those regions that contain e and take a vector £ starting at vq and 
ending somewhere on the intersection 

n n «j) + 

v a vertex of P j s.t. 

(e—v,a v j)>0 

(we can take for instance £ = e — vo); then clearly (£,a v j) > 0 for all edge vectors 
a v j satisfying (e — v,a v j) > 0 (see Figure|Gl). □ 



Figure 6. 

Remark 38. We have seen in the above proof that, choosing e in an appropriate region, 
the polytope decomposition formula I® becomes the Lawrence-Varchenko relation. 
In addition, in some cases, we can also choose e so that El becomes the weighted 

^This weighted version of the Lawrence-Varchenko relation is different from the ones in IKSWI 
and EU because here we may assign different weights qi E C to the faces of the polytope instead of 
a fixed complex number. Nevertheless, the proof of this decomposition formula follows easily from 
the ones in IKSWI and ixn by applying Lemma Idol to the boundary points, as we do in the proof of 
Theorem rra 
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Brianchon-Gram formula of IA2I . Indeed, considering for each vertex v of P, the cone 
generated by the inward normal vectors to the facets through v, and taking the 
intersection 

p d ■■= n c - 

v vertex of P 

then, whenever int(P <2 HP) ^ 0, we can take e G int(Pd HP), and obtain toa = 0 and 
tp(e, A) = 1 for every A in B. Then, with this choice of e, (13011 becomes the weighted 
Brianchon-Gram formula: 

(39) lp = E(-l) dimF lc F > 

F 

where the sum is over all faces F of P. 


5. The weighted Euler-Maclaurin formula 


As an application of our polytope decompositions, we will give new weighted Euler- 
Maclaurin formulas with remainder for the sum of the values of a smooth function / 
on the integral points of a simple polytope P. 

5.1. Weighted Euler-Maclaurin for intervals. Let us first recall the weighted 
Euler-Maclaurin formula for this sum presented in |AW| (see also |Kn| and jKSWj l: 
let q be any complex number and let / be any C m function on the real line (m > 1); 
for integers a < b, the sum 

(40) E 9 / := «/(«) +/(“ + !) + ••• + f(b-l)+qf(b) 

[a,b] 


rb+h.2 


= Qf(Pi)Qf(P 2 ) 


/( x) dx 


' a—hi 


h±=h2 —0 




where k = |_m/2j, where Q q k (S) denotes the truncation at the even integer 2k of the 

power series 

(41) 


Q,(5) = (g-l)5+Td(5) = 1+ ( g - E P+E 77^7 S 2k =(q- 


s+ 


S/2 


tanh(5/2) 


(here Td is the classical Todd function defined by Td(S) := S/ (1 — e~ s ) = 1 — biS + 
SfcLi jWj\ ^ 2fc ’ bk the fc-th Bernoulli number 0]), where 


and 

(42) 
with 

(43) 


n d _ d 
1 dV 2 dh 2 ' 


Rm(f) ■= (-l)™" 1 [h P m (x)fW(x)dx, 

J a 


P 2k+1 := (-I)*’ 1 E 


2 sin(2?r7rcc) 


^ (2nn) 2k+1 (2k + 1)! 


P 2 fc+l({^}) 
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if m = 2k + 1 is odd, and 


(44) 


P2k 


(-i )^E 

n= 1 


2 cos(2mnc) 
(2n7r) 2fc 


1 

m 


B 2 k({x}) 


if m = 2k is even, (here B m is the m-th Bernoulli polynomial and {a;} := x — |_£j is 
the fractional part of x ). 


Remark 45. The functions satisfy the following symmetry property 

(46) Q ?{S) = Qi- q (~S). 

Indeed, Qq k (S) is a polynomial with constant coefficients, 1 + (g — ^)S+ terms of 
even degree independent of q. 


Equation (in . when applied to a C m function of compact support, gives the 
weighted Euler-Maclaurin formula for the half ray [a, oo): 

(47) 

/•OO 

^2 q f '■= if (a) + /(« + 1) + f(a + 2) 4 -= Qjf (£>i) / f(x) dx 

\ J a—hi 


hi —0 




[a,oo) 

where 

/•OO 

(48) R m := (-I)™" 1 / P m (x)f^ m \x)dx. 

Moreover, for the half ray (—oo, a], we have 

(49) q f = qf(a) + f(a - 1) +/(a - 2) + ••• 

(—oo,a] 

and so, considering the function g defined by g(x) = f(a — x), we obtain 

X) 

g(x) dx 


(50) '52 q f='52 q g = Ql k (D 1 

( — oo,a] [0, oo) 

r-a-\-hi 


I-hi 


= Qq ( D l) 


/( x) dx 


hi=0 


poo 

+ (-l ) m - 1 Pm(x)g {m) (x)dx 

u=o J o 

+ Rm(f), 


where 

(51) 


Rm(f) := (-I)" 1 ’ 1 f P m (x)fW(x)dx 

J — OO 


(here we used the parity and the 27r-periodicity of sin(x) and cos(a;)). From formulas 
m and m and symmetry property (gnji, we obtain the Euler-Maclaurin formula 
for the whole real line K: 
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5.2. Twisted weighted Euler-Maclaurin formulas for intervals. We will now 

consider the twisted weighted sum for a half ray 

OO 

(53) E 9A "/w = 9 /(o) + E A "/( ri ) 

n>0 n=1 

where A ^ 1 is a K-th root of unity with K a positive integer. Let Qm,\ be the 
distributions defined recursively in IKSWI by 

Qo,\(x) -=~Y An ^( x _ n ) 

nG^ 


and 


7 pK 

—j — Qm— i,a(*^) and j (Lx — 0. 

ax J Q 

Moreover, let us consider the polynomials defined in AW by 

N k q ' x (S) := (q + S + Q 2 ,a(0)S 2 + Q 3 ,a( 0)S 3 + • • • + Q fc ,A(0)S* 

where A ^ 1 is a root of unity. 

Since ^lr„ iTl _|_i)(x) = 5{x — n) — S(x — (n + 1)), we have 

E (^A"i Kn+1) (x)j =EdEA^(*-«) = ^Q 0l A(*), 

VnG^ J 


K -1 


implying that Qi,\(x) = Yj n ez ^ ni [n,n+ 1 ) (note that f Q Qi,\(x) dx = T ^ ^ A’ ; 

71=0 

0). On the other hand, integrating by parts, we have 

poo \ oo pn-\-l \ 

J o Qi,\(x)f(x)dx = Y^Ey X n f'(x)dx = - —^/(0)+A/(l)+A 2 /(2)+- 
and so, 

A r 00 

g/(°) + E A ”/( n ) = (9 + Y3^)/(°) + J Q QiA x )f'( x ) dx 


n> 1 


= (g + E)/(0) - Q 2 , A (0)/'(0) + q 3 ,a( 0)/"(0) - • • • + (-l) fe - 1 Q fe ,A(0)/ (fe - 1) (0)+ 

poo 

+ (-l) fe_1 / QkAx)f {k \x)dx. 

Jo 

Then, since (— l) m f( rn ~ 1 ' > (0) = (^) m fA f( x ) dx\ h=0 , we obtain the following twisted 
Euler-Maclaurin formula: 

Proposition 54. (!Mj, IKSWI ) Let k > 1 and let f £ C fe (KL) be compactly supported. 
Then 

(55) Y , ,A "/(n) = N*’ A f JA r f{x) dx + r Q M (*)/(*>(x) dx. 

Eo /,= 0 Vo 
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Remark 56. If, for A ^ 1, we write A = e 2m ^ K , then, by the Poisson formula, we 
have 


Qo,a(x) = — ^ A” 8{x — n) = — ^ A x <5(:r — n) 


n£Z 
„2nix-z? 


n£Z 


£ ^ 6(x - n) = -e 2 ™7 Y, e ^ lTX = ~ E 1 


,27ri(r+-£.) x 


nEZ 


re Z 


rez 


Hence, for ?n > 1, we obtain 
(57) 


1 e 27ri(r,+ ^ )a: 
Qm,A\X) TZ TT~” / ' T~ 5 

( 27r *) ^ (r + je) m 


and so Q m ,A(0) = — Z)rez ( r+ j ^ m is the (m — l)-th coefficient of the Taylor 


series expansion of 




= Y 3 J- — at s = 0 (the derivative of 


respect to s is equal to 


= i r 

4tH 


l_e 27 ri i-- s 


with 


-, and higher order deriva- 


4 sill 2 (2J--4L) *** L-irel ( r + i _^)2 

tives are obtained differentiating this series expansion). Consequently, considering the 
operators 

T(A ’ s): =d7 

defined in IBrVlI . we have that N q' X (S) is the truncation at the integer k of the power 
series 

N 9 A (S) := (q + Y ^j)S - YZT\ + T ( A ’ S ) = (v- !) 5 + T ( A > S). 

From eo> it is clear that the operators N™ ,A satisfy the following symmetry prop¬ 
erty 

(58) 


\s) = N™’ a (-S'). 


Remark 59. If, for A = 1, we define 

N k /(S) :=Q 2 L fe / 2 J(S) and Q k ,i := -Pfc, 

then formula <ESJ» becomes formula fn and so it is still valid. Note that, if A ^ 1, 
N a ' ,a ( 5') is a multiple of S and that, if A = 1, then N^ ,a ( 5') = 1+ a multiple of S. 
Moreover, still when A = 1, symmetry property becomes property {El)- 


5.3. Weighted Euler-Maclaurin formulas for cones. For a subset Jc {1,..., d}, 

let S j be the standard J-sector S j := {x G l d Xj > 0 for j G J}. Iterating 
equations 63 and iTrl . we obtain an Euler-Maclaurin formula for S j (J ^ 0) and 
a C m function of compact support: 


(60) 

E 

Sj nz‘ 


7:=E 

Xj EZ+jEJ 
Xj EZ, j(£j 


/it..«i = n^) 


'Sj(M 


/(x) dx 


hj= 0 


■^m(/), 


^Note that -r-£- 

sin-^ 


E/rEZ (r+z)^ 
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where lg is the weighted characteristic function for the J-sector defined in Re¬ 
mark Em where Di = d/dhi , where hj = {hj 1 ,..., hj n ) with J = {j i,..., j n }, where 
Sj(hj) = {a: G | Xj > —hj , for j G J} is the shifted J-sector, and where the 
remainder R^ t (f) is given by 


RtU) ■■= 


= EE( _1)(to-1)(|-R 


RC{l,...,d} 

R=£I 


If J = 0 then Sj is the whole space and so 


I —171 


IlQ 

iei 


«(A) 


/ ip^ii 

'Sj(hj) ieR \j jeR\f 


dxj 


f(x) dx 


hj=0 


(61) 


E w / = 

Sj nz d 


f(x)dx + R®(f), 


with 

(62) <(/) 

Let us now consider a regular integral J-sector Cj, the image of the standard J-sector 
by an affine transformation 


E 

R^0 

RG{l,...,d} 

R^I 


(-1) C 


m — 1) |i?| 


i£R 


j&R 


n p m {xi) n ( dx 


x i—> AjX := Mx + b, with M G SL(d , Z) and b G R d . 


Moreover, let us denote by C j(h) the expanded sector, image of S j(h) under this 
affine transformation. For a C m function of compact support /, let us consider g := 
A *jf = f °Aj. Then, 


E 

Cj nz 


’/:= E 

Sj nz' 


3=IlQ 

j&J 




ISj(hj) 


g(x) dx 


hj=0 


+ Rfr (s), 


and we obtain the following Euler-Maclaurin formula for a regular J-sector: 


(63) 


\ ’ 


/=nw/ f^ dx 


c j n i 


jeJ 


'Cj(hj) 


hj —0 




where R° J (/) = R^{g). 


5.4. Weighted Euler-Maclaurin formula for regular simple integral poly¬ 
topes. From © we can write an Euler-Maclaurin formula for a regular integral 





NEW POLYTOPE DECOMPOSITIONS AND EULER-MACLAURIN FORMULAS 


23 


polytope P with N facets, by using a polytope decomposition from Theorem 14.11 


E"7:= E 1 p/=E(- 1 ) m M £ > A ) E “/ 

pnz d pnz d AgB cj^nz 1 * 



where 

(65) S?(f) := E(-l) m Ma )R${f), 

Agb 


and where the dilated polytope P(hi ,..., /ijv) is obtained by shifting the itli facet out¬ 
ward by a “distance” hi. Here we used the fact that, when multiplying the differential 
operator in the first term of the right hand side of (p|) by any operator of the form 
Q 2 q k j(Dj)i with j (ji J, all that will remain of Q q k {Dj) is the constant term 1, not af¬ 
fecting the final result. Note that both ^ pn z” / and IlyLi fpfa h N ) f( x ) 

do not depend on the choice of e (that is, do not depend on the Paradan region used). 
Consequently, the remainder is also independent of this choice. 

Remark 66. Alternatively, using the polytope decomposition of Theorem l4.2l we obtain 
a different expression for the remainder in (|65l). Indeed, we get 

(67) $£(/) := E(- i r A+dimA ^,A) J R^(/), 

AgB 

where now the tangent cones Ca are polarized with respect to the vectors e — j3(e, A). 

5.5. Weighted Euler-Maclaurin formula for simple integral polytopes. To 

extend formula (IMll to simple integral polytopes we need to obtain an Euler-Maclaurin 
formula for simple J-sectors. We can describe a simple J-sector, Cj, with J = 
{*i,..., i„.} C {1,..., d} as the intersection of n half spaces Hi in general position 

(68) C j=f] Hj , 

jeJ 

where Hj := {x £ | (x, ijj) + A j > 0} for rational vectors rjj. Clearly, C j is a cone 

of the form along the affine space A j S B defined by 

At = p| dH J 

jeJ 

(see C5» . Clearing denominators we can assume the r/,; ! s to be integral and we 
impose the normalizing condition that they are primitive elements of the dual lattice 
Z™* (note that these vectors are inward normals to the facets of C j). Let us take 
the dual basis {aq,..., ai n } in R” (that is, such that (a*,, rji) = Ski for k,l € J) and 




24 


J. AGAPITO AND L. GODINHO 


denote by Tj C the subspace generated by these vectors. The oti s are what, in 
Section^ we called generators of Caj and generate a lattice i in K"* which is a finite 
extension of Z ra (this extension is trivial exactly when Cj is regular). Let Tj be the 
finite group 

Tj = (Z"* n Tj)/t. 

This group is trivial exactly when C j is regular, and its order is |Tj| = | det t\. 

Moreover, as it is shown in IKBW i. 7 1 —> defines a character of Tj, whenever 

x £ £, which is trivial iff x £ Z ra . Since, by a theorem of Frobenius, the average value 
of a character on a finite group is equal to zero if the character is non-trivial, and 
equal to one otherwise, we have 

1 y- e 27ri<7,x) _ ( 1 if x £ Z" 

(TjI A I 0 otherwise 

7 fci J V 

for all x £ t. Consequently, for any compactly supported function / on R n , 

(69) £“/ = 

c jnz d ' 7er x 

where we sum over all 


(70) 


= y + J2‘ 




with y £ A ,7 D Z ra and all rrij £ Z + with j £ J. Moreover, the cone C j is the image 
of the standard J-sector Sj under an affine map 

(71) 1Ajt := Ujt + b, with b £ l d , 

where, for j £ J , Uj £ GL{d 1 Z) carries the vectors ej of the standard basis of W 1 
into the basis {a^}]J =1 . Hence, \detUj\ = 1 /1 det= l/|Tj|. On the other hand, 
since in cm we have y £ Z”, we get 


0 2iri(y,x) _ 


ro; 


with Xj = e 27ri( ' J ’ aj ), 




and so the inner sum in m becomes 


E E E* n^' /fe+E= 

jEJ yGAjflZ 71 rrij >0 \j£j / 


ieJ 


(72) 




j£J l<i<d mi^L rrij >0 \lGJ 

i£J 


where g = f o Aj. Iterating the twisted remainder formula for the half ray 
and the Euler-Maclaurin formula for the whole real line, the sum in m can be 
written as 

d 


(73) 


n 


dh, 


f Sj(hj 1 ) 


gj(t)dt 


h —0 


+ -Rq r,fcC\jl> ■ ■ • ) ^jn’l9 )) 
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where again S j(hj 1 ,... ,hj n ) = {(ti,...,t^) | ti > —hi for j £ J} denotes the 
dilated standard J-sector, and where, for qj := {qj 1 , ..., qj n ), the remainder is given 

by 


(74) < d fc (A 7 j 1 ,...,A 7ljB ; ff ) := £ £ 


ICJ KDJ 
RC{l,...,d} 
R=£I 


n 


iei 


dh, 


ij -/s/;,. h 


IT (tj) ]^[ 


jeR\! 


jGR\I 


dt 1 


git) dt 


h =0 


with g = f o A. Changing variables by the inverse transformation of CD- the Euler- 
Maclaurin formula in (H3> becomes 


(75) 


c J nz d 7 er ieJ 



'Cj(hj) 


f(x) dx 


hj=0 


+ <%(/) 


where A 7J := e 2 ’ r ®<7,aj) ) where, for hj := ihj lt ..., hj n ), Cj{hj) denotes the image of 
the dilated standard J-sector Sj(/ij) under the affine transformation Aj defined in 
CD. and where the remainder fc (/) is given by 


(76) 


R Cj 

U ci j,k 


if ) : =EE E (-i) (fc - l)(|fl| - |j| )n^-‘ 

7 er/CJ RDJ iei ' 1 

R=£I 



b )) n D jfi x ) dx 

jeR\i 


h—0 


with iU-^k the jth row of U 1 and with Dj the directional derivative along the jth 
column vector of £/ -1 . Note that, when j £ J, this is the directional derivative along 
a j- 

Let now P be a simple polytope and choose an e on some Paradan region. For each 
affine space A generated by a face of P there is a J-sector C j equal to the polarized 
tangent cone of A (cf. itTTI) ') and so we can associate a finite group Ta to A 
by simply taking the corresponding group Tj. Let P{h ) denote the dilated polytope 
obtained by shifting the *th facet by a distance hi. Our decompositions of Pih) involve 
dilated sectors but now, dilating the facets of P outward results in dilating some of 
the facets of inward and some outward. Explicitly, taking J C {1,..., d} such 
that = C j (see CHJ)), the inward normal vector to the jth facet of (j £ J) is 

J j Vj, if <*aj = 

Aj 1 -r tt > 

[ ~Vj, “A ,j = -«A ,j- 

where qj is the inward pointing primitive normal vector to the jth facet of P (note 
that <Ta j G J, is the dual basis to the corresponding vectors qj). The dilated 
sectors that appear on the right side of the polytope decompositions of P(h) are then 
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C A( h A ■ • • j h A,j n ) ( with J = ij i> • ■ ■ ’jn}), where 


fot _ I ^ a A,ji — a A,ji 

A ’ Ji \ ~ h 3ii if a A,ji = '“A ,n- 

Moreover, the roots of unity that appear in the Euler-Maclaurin formula for are 

\tt _ „27Tj<7 ,ott ,) _ J A ' if a Aj = a A,j 
A 7 j,A — e ,J — S X -1 .j. )j 

l A 7 ,j,A> lf «Aj = ““A ,j. 

Hence, for any compactly supported function in R d of type C nk (for an integer 
k > 1) the decomposition formula of Theorem 14. 1 1 annlied to P(h) along with formula 
(1751) give 

( 77 ) E w f= E(-irx^A) E w f = 

pnz d AeB cl nz" 


E(-irM^A)E n n 

a eB 7er A jeJ A 


k, A" 


d 


7 _ 

?! l 


Aj7 -M(^) 


/(x) dx 


ft," =o 


+ Rw,k(f)> 


where h j = (ftj^,..., h\ n ) and where 


(78) 


<fc(/) := E(- 1 ) mA ^ e > A ) i? 5 A 


A<E£ 


Remark 79. Using the polytope decomposition of Theorem 14.21 we obtain 

(so) E “/= E(- 1 ) mA+dimA ^ e ’ A ) E ”/ = 


pnz d 

E(_ 1 ) "*A+di m A (£)A) E n N fcAUA f-§- N ) / t f i x ) dx + <*(/), 

Ja={ji, —,jn A } 

with h'j = (/i^,..., /i]j ) and 

(81) <*(/) “ E(- 1 ) mA+dimA ^(^ A ) J R? ,(/), 

AeB 

where now the tangent cones Ca are polarized by the vectors e — /3(e, A). 


Aes 


c*nz" 


Let us now analyze some properties of the groups Ta- These generalize Claims 61, 
62 and 65 in jKSWj to spaces A G B of arbitrary dimensions. Their proofs follow 
easily from the ones in IKSW| but we will include them for completeness. For that 
we will first introduce some necessary notation. If A and A are two elements of B 
with AC A, then clearly C Ta- Hence we can define a subset T^ of Ta by 

r i :=r A \ (J r A 

AeBj ACA 
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and then 
(82) 


r A = U r A 

AeB| ACA 

Claim 83. If j £ Ta and j £ Ja, then Xjj.a 1 is the same for all A' C A. 

Claim 84. If 7 € Ta, A' c A and j £ Ja' \ Ja, then A 7 J! a' = 1. 

Claim 85. If 7 £ r A and j £ Ta, then A 7J ,a 7 ^ 1. 

Proof. Let 7 £ Ta be represented by 


7 = X! biT]i G Tj a n 

iGJa 


c£\* 


for some bi £ R. Let A' C A. Since the oia'j’s are dual to the rjj’s, for j £ Ja>, we 
have 

bj, if j £ Ja 
0 , if j £ Ja’ \ Ja 


(7, OtA-,j) = 


Consequently, 


A, 


■•yJA' 


e 2 * ib i , if j £ J A 


1, if j £ Ja' \ Ja 
is independent of A' and is equal to 1 if j £ Ja 1 \ Ja, and so Claim OSH and 1841 follow. 
Let j £ Ja- If A 7 , ?i a := e 2mbj = 1 , then bj £ Z and so 

( 86 ) 


7 = ^ 


also represents 7 . Let A D A be the element of B such that J A = Ja \ {j}- Then, by 
1861) . 7 £ r A , and Claim 1851 follows. □ 

With these properties we can further simplify formula O- First, note that either 
h A,j = h P A 7 .j,A = A 7 ,j,a and q) = q 3 , or h} A j = -hj, A^- A = \~] A and q) = l-q 3 , 
and so, by symmetry property (1551) . this gives 


(87) 


fc,A # a ( d 

N 7 ’ 3,A * - 

< j 




Hj 


_d_ 

dhj 


Moreover, from Claim l%H we have A 7Ji a = 1 for j (f Ja, implying that 


K ^ { JL ) = 1 + powmo! A. 


Since, still for j ^ Ja, the cone C A (/i A ) is independent of hj, ( 1771 ) is equal to 


N 


(88) ^(-irMc,A) £ II IV A ; A - - 


A 60 


tsTa j=i 


_d_ 

dh. 


■iJ 


f(x) dx 


h =0 


+ Rw.k(f), 
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where N is the number of facets of P. Defining 

N / d \ 

(89) N* A : =n N £ W ( Wj )> for 
we have, from Claim IH3 that 

(90) N^ A = N fc ~ whenever 7 G T A and A C A. 
Consequently, using 15^1 . formula dSHl can be written as 


(91) 


£(-irMe,A) E N 7 ,a , f(x)dx +<*(/) 
AeB 7 er A h =° 


= E(- 1 ) m M £ ’ A )E EE 


/( x) dx 


AeB 


AeB 7ert 




+ KM) 


= E E^iEHrx^)/, ,/(*)<& +<*(/)■ 


AeB 7 er 


AcA 


In the interior summation on the left we can add similar terms that correspond to 
spaces A not included in A. indeed, these make a zero contribution for the following 
reason: if A is not a subset of A then there exists a j G J A \ Ja; then since j £ Ja, 
the cone C A (/i/ A ) does not depend on h } ; on the other hand, since 7 G and 
j G J A , we know, from Claim[S51 that A . A ^ 1 and so, by Remark 1531 we have that 
N q j Xy ' :l ’ A (one of the factors of ~) is a multiple of (^-). Therefore, (ITJTll is 

equal to 


E EEiEHrM^) 

AeB 7 er b AeB 

(92) 

= E E n E 

AeB 7ert 

and we have our result: 


!p(h) 


f(x) dx 


cUh*j) 


/( x) dx 


h=0 


h =0 


+ K,M)’ 


KM) 


Theorem 5.1. Let P be a simple polytope in R d with N facets and let f G C)? fe (K d ) 
be a compactly supported function on for k > 1. Choosing an e on a Paradan 
region determined by P, we obtain 


(93) 


E*/=E E N C 

pnz d AeB 76 r b i 


lp{h) 


f (x) dx 


h —0 


K,M), 


where A is the differential operator described in m and TO . and where the 
remainder is given by ca. The operator N* A is of order < k in each of the variables 
hi,... ,/ijv with N the number of facets of P. The remainder is a sum of integrals 
over sectors, of bounded periodic functions times several partial derivatives of f of 
order no less than k and no more than kd. Moreover, this remainder is independent 
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of the choice of Paradan region of e, and is a distribution supported on the polytope 

P. 


Remark 94. If we instead use the polytope decompositions of Theorem 14.21 the re¬ 
mainder in El will be given by ED. 


The Euler Maclaurin formula El obtained in Theorem rm is similar to the one 
presented in rawi . However, in our formula, we allow the operators Ngj 7 ’*’ A that 
define A and Nq to have different weights qj £ C while, in rawi . the qj's are all 
equal to some fixed complex number (in IKSVVI this fixed weight is 1/2). Moreover, 
we obtain a different expression for the remainder k (f) which is now given as a 
sum over the affine spaces generated by all the faces of the polytope (not only over the 
vertices). In addition, the intermediate formulas that we obtain in (191 1) (before adding 
terms with zero contribution in order to get an integral over the dilated polytope) 
also involve sums of integrals over the polarized tangent cones to the polytope at the 
different faces and not only at vertices. 

Just as the Euler Maclaurin formulas in rawi and KS Wl our formulas generalize 
to symbols 5 , giving rise, in particular, to the following exact formula for a polynomial 
function p in K. d 


(95) 


E"> 

pnz d 


E E N U 

AeB 7er) i 



p(x) dx 


h —0 


(where we choose k > deg p + d + 1). From Remark ESI we see that this formula is a 
weighted version of the exact Euler Maclaurin formula obtained in IBrVlI . which is 
obtained from ED by making all the weights in w equal to 1. 

Acknowledgments: We would like to thank Jonathan Weitsman for calling our 
attention to the localization method of Paradan applied to the norm square of the 
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